DYNAMIC EQUIVALENCE OF CONTROL SYSTEMS VIA 
INFINITE PROLONGATION 



MATTHEW W. STACKPOLE 

Abstract. In this paper, we put issue of dynamic equivalence of control 
systems in the context of puUbacks of coframings on infinite jet bundles over 
the state manifolds. While much attention has been given to differentially 
C flat systems, i.e. systems dynamically equivalent to linear control systems, 

^ the advantage of this approach is that it allowed us to consider control affine 

systems as well. Through this context we are able to classify all control affine 
systems of three states and two controls under dynamic equivalence of the type 
CN (x, u) i-s> y(x,u). 



1. Introduction 



q 

'"pi A control system is an underdetermined system of n ordinary differential equa- 

C3 tions (ODEs), 

, (l.I) X = f(x,u). 

Control systems show up in the design of electrical and mechanical systems, among 
" ^ other things. The variables x whose time evolution is determined by the ODEs are 

— called state variables, while the "free parameters" u are called control variables. A 

control system can be viewed as a submanifold S of the tangent bundle of the state 
space in the following way: given a manifold M and a curve x : / — )■ M, we say 
^ that X is a solution to the system S C TM if (x(t),x(t)) hes in E for all tel. The 

\^ map — > TxAf given by u H> (x, f (x, u)) is a parametrization of Sx = E n T^AI 

with the parameters u seen as local coordinates on Ex. 
] A dynamic equivalence takes trajectories of one system, x = f(x, u), to those 

. . of another, y — g(y, v), and back again via maps between jet spaces which allow 

^ state derivatives to get mixed in: 

X (x, u, u, . . . , u^'-')) y(x, u, u, . . . , u^-')). 

Through the defining equation , derivatives of state variables can be expressed 
in terms of control variables and their derivatives as well. Static (feedback) equiva- 
lence, which is a diffeomorphism of the state space, is a special case where y = y(x). 

Up to dynamic equivalence at the first jet level (J = 0), i.e. x — x(y,v) and 
y = y(x, u), my results classify all control affine systems, 

x = f°(x) + u,r(x), 

of at most three states and two controls through the use of Cartan's method of 
equivalence. The main result of this paper is that every control affine system of three 
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states and two controls falls into one of three classes under dynamic equivalence. 
The numbered rows represent these three classes. The entries in each row are 
systems that, while dynamically equivalent, are not statically equivalent. 



1 


Xl — Ul 


Xl 


= Ul 


Xl — Ul 




X2 — U2 


±2 


= U2 


±2 = U2 




is = X2 


±3 = 


- X2U1 


±3 = 1 + X2U1 


2 


Xl = Ul 
X2 = U2 

is = 








3 


Xl = Ul 
±2 = U2 
X3 = l 









2. Control Systems 

A system of ordinary differential equations (ODEs) with more variables than 
equations is called a control system. Locally a control system with n+s variables 
and n equations can be written in the form 

Xl = fl{xi, Xn, Ul, Us), 

X2 = f2{xi, Xn, Ul, Us), 

Xn ~ fn{xi, ' • • , Xn, Ul, • • • , Us)- 



For our purposes, we will consider the functions /j, 1 < z < n, to be C°°. 

Here, : M — >■ M and uj : M. ^ M. We will use t as our independent variable, 
and derivatives with respect to t will be denoted by a dot: ^ = .Xj. This system 
of equations can be abbreviated with the single vector equation x = f (x, u) where 
X = {xi, Xn)'^, u = {ui, Us)"^, and f = (/i, fn)"^- This type of 
control system is called time independent since there is no explicit t dependence 
in the /». 

In general, quantities that are vectors or matrices, like x above, will be written 
in bold face to distinguish them from scalars. like .t,. 

The variables known as the state variables, while the variables uj are 

known as the control variables. To explain the terminology, imagine a hovercraft 
on the surface of a two-dimensional lake. The state variables would be those needed 
to describe the state of the hovercraft on the lake: position of the hovercraft, which 
direction the hovercraft is turned, and the translational and rotational velocities 
of the hovercraft. The time evolution of state variables is predetermined, in this 
case by the Newton-Euler equations of motion, which are given explicitly in the 
example below. The control variables allow external influence of the state variables' 
time evolution. In the hovercraft scenario, control variables could describe the 
hovercraft's motor: the magnitude and direction of its thrust. Control variables 
are exactly what the hovercraft operator uses to control the system. 
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3. Dynamic Equivalence 

Geometrically, a control system can be viewed as a submanifold S = ]R x S of 
M X TM in the following manner: Given local coordinates x on M , the control 
system S is a manifold with local coordinates (x, u). With local coordinates (x, x) 
on TM, there is an embedding 

i, : ffi X S K X TM 

given in coordinates by 

(i,x,u) ^ (t,x,f(x,u)). 

This embedding /, pulls back the contact forms {dxi — Xi dt \ i = l,...,n } on 
R X TM to the forms { uj, = dx, - /,(x, u) | i = 1, . . . , n } on M x S. 

A solution to a control system also has a geometric interpretation. Let x(t) be 
a curve in M, i.e. x : M M, and define pix(t) = (x(i),x(t)) e TM. Such a 
CTirvc x(t) is a solution to the control system E if there exists a map cr : M — ^ S 
that makes the following diagram commute: 



IS 

E ^TM 




M 5>M 



In particular, 

Pix{t) = o cr) (t), 

or in other words, pix{t) G E for all t. Note that uJi{&) = for i = 1, . . . , n. 

We will use the convention in this paper that a control system with a bar over 
it, for example E, is a subbundle of K x TM, while a control system without the 
bar, E, is a subbundle of TM which is the projection of E. In fact, since we will be 
requiring that time be preserved through our equivalences, we will have E = M x E. 

3.1. Jet Spaces. Since the idea of dynamic equivalence is to allow a "change of 

variables" using higher order derivatives, wc need a setting in which these higher 
order derivatives can be dealt with, much like the tangent bundle lets us work 
with first order derivatives. This setting is a jet space. We will say that curves 
a, 6 : M K with a(0) = 6(0) = have the same K-jets at if 

da db d'^a d^b d^a d^b 

di^'^=dt^'^' ^W = ^(0)'-'^(°) = ^(0)- 
Given n-dimensional differentiable manifolds U and V and maps a,b : U ^ V with 

a{x) = b{x) = q, we will say that a and b have the same K-]ets at x if for any 
differentiable maps (p : M. ^ U, ip : V ^ M. with (j){0) = x, tp o a o (j) and ijj o b o cp 
have the same JT-jets at 0. 

Note that having the same K-}ets at x is an equivalence relation among maps 
from U to V. Define the ii:*'*-order jet bundle of M, denoted by J^{M), to be 
the bundle over M whose fiber i7^(M)x over a point x e M is the space of curves 
a :M. ^ M modulo the equivalence relation of having the same ii'-jets at x. Notice 
that with this definition, J°{M) = M and J^{M) = TM, where the equality here 
is actually a bundle-preserving diffeomorphism. 
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Define the prolongation map pj^k wliicli takes lifts of C°° curves from M in 
Ji{M) to lifts of C°° curves from M in j''{M) {j < k) as follows. 

p,- ,(x(t),x(t),x(i), . . . ,xW(t)) = (x(<),x(t),k(t), . . . ,xW(t), . . . ,xW(t) ) 
We will denote pQ j simply as pj. 

3.2. Definition of Dynamic Equivalence. Let M and N be smooth manifolds 
(state spaces) and 



(3.1) 



E : X = f(x,u) 
A : y = g(y,v) 



control systems over their respective state spaces. 



We say control systems (3.11 on M and N are dynamically equivalent over 



open sets U C J''^^^{M) and V C J^^^^^{N) for nonnegative integers J and K if 
there exist smooth maps ^ : U ^ N and 'J : V -> Af so that when restricted to 
the appropriate open sets: 

(1) for any solution x(i) of x = f(x, u), ($ o pj^i)[x{t)) is a solution to y = 

g(y,v), 

(2) for any solution y(<) of y = g(y,v), (^f o PK+i){y{t)) is a solution to 
X = f (x, u). 

(3) the following diagram commutes for solutions, 



PK+iy 




i.e. '5 o px^i o $ opj^i{x{t)) = x{t) whenever x{t) is a solution of E, 
and $ o o v]/ o j)/f_|_i(y(i)) = y(i) whenever y(i) is a solution of A. 

Note that this means 



= $(^X,i,.. 

= *fy,y,- 



,y 



We will use the same notation for maps between jet spaces as we did for control 
systems, namely ip : J^{M) J^{N) and : M x J^{M) ^ M x j''{N) with 
If = id X if. Also note that in the definition of dynamic equivalence, we are using 
maps if : J'{M) J^{N), so they are defined in terms of the coordinates: 



( X, X, X, 



x(^) ) ^ ( y, y, y, . . . , y^^'^ ). 



(fc) 



However, in practice we will be concerned only with the restrictions of these maps 
to the prolongations of control systems (defined below). Therefore, by way of the 
defining equations x = f(x, u) and y = g(y,v) of the control systems, we will be 
looking at the restriction of Lp to the appropriate submanifolds with the following 
coordinates: 



( X, u, li, 



u 



(j-i) 



) 1^ ( y, V, V, 
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The proof of the following theorem should be clear from the definition, which is 
the same definition given in |10| . 

Theorem 1. Dynamic equivalence is an equivalence relation of control systems. 

Static (feedback) equivalence is a special case of dynamic equivalence for 
which J = K = —1, i.e. $ : Af — > is a diffeomorphism with = For static 
equivalence, we have 3>*E = A and — S. We say two systems are locally 
static equivalent over U C M and V C A^ if there exist coverings t{ = UasA^" 
and V = Uqsa such that the systems are static equivalent over each and 
V". 

From an engineering point of view, equivalence can be achieved through the 
addition of a feedback loop in the control system. In Figure [T] the system E has 
input u and output x. By adding a feedback loop, the new system A has input v 
and output y. In the case of static equivalence, the feedback loop only incorporates 
the old output x so that the new input v is only a function of x and u. Including 
one or more integrators to the feedback loop allows v to be a function of x, u, and 
some number of derivatives of u, and this is dynamic extension. 



A 



u 



Figure 1 . Control system with feedback 



3.3. Prolongation. A key ingredient in dynamic equivalence is the notion of pro- 
longation of a control system. For integers fc > 1, define the prolongation of 

the system S to the fc*'* order, denoted by S^, to be the subbundle of J'^{M) that 
corresponds to the prolongations of solutions of E, i.e. for any x : I ^ AI , 

Pi(x(t)) e E Vt e / ^ pfc(x(<)) e Efc Vt e /. 

Obviously Ei = E. In the same way that E is a control system with s control 
variables with state manifold M of dimension n, we can view E2 as a control system 
with s control variables with state manifold E of dimension n + s. An important 
fact is that E is strictly dynamically equivalent, i.e. dynamically equivalent 
but not static equivalent, to E2, as can be seen in the diagram below. 



6 



MATTHEW W. STACKPOLE 




J\M) 

Example 1. The system 

Yi : xi = ui 

has two states and two controls. E is dynamically equivalent to 'E2- 



yi 
m 

2/3 



Vl 
V2 



where 
(3.2) 



Xi 
U2 



yi, 



X2 

ill 



y2, ui 

Vl, U2 



2/3, 
V2- 



We have increased the number of states from two to four by viewing the controls as 



new state variables. (3.2) gives the equivalence map. This is an example of what 



we will call a total prolongation. 

In general, a total prolongation of the system x = f{x, u) is the system 
X \ / f(x, u) \ , . 



(3.3) 



u 







where is the vector with a 1 in the {i + n)*^ entry and zeros elsewhere. Here 
(x, u) are the new state variables and u are the new control variables. This system 
has a special form. A control system of the form 



(3.4) 



x==f(x,u) =f°(x)- 



r(x)u. 



is called control affine. In particular, (3.3 1 is control affine. Thus we have the 
following theorem. 

Theorem 2. Every control system E is dynamically equivalent to a control affine 
system, namely E2. 

Similar to a total prolongation, some, but not all, of the control variables can be 
made into new state variables, as we see in the following example. 
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Example 2. The system 



is dynamically equivalent to A: 



where 



Xi 
U2 



yi, 



E : xi 
A: yi 

m 



2/2, 
V2- 



Ul 

2/3 
Vl 
V2 

Ul = 2/3, 



We have increased the number of states from two to three by viewing only one of the 
controls as a new state variable. This process is called a partial prolongation. 
Every control system is dynamically equivalent to any partial prolongation of that 
system. 

We will assume, without loss of generality, that the two systems in a dynamic 
equivalence have the same number of state variables (to = n). If to < n, perform 
repeated prolongations, either partial or total, until the number of states are equal 
and consider this new system. 

A method for constructing a potential dynamic equivalence which is not a partial 
prolongation was mentioned briefly in a paper by Pomet Below we give a 
specific example of how the method works. This example incorporates both partial 
prolongation and changes of variables (a.k.a static equivalences) to give not only two 
control systems that are strictly dynamically equivalent but also the equivalence 
map. 

Example 3. Start with an affine linear control system: 





(3.5) 



Xl \ 




( 1 


±2 


- 





±3 / 




I X2 




U2- 



Partially prolong the three state system to a four state system. 



Zi — Xl 

zs. = 

Wi = Ul 



Z2 = X2 
Za = U2 
W2 = U2 







1 






( ' ] 






Z2 


















Z3 









+ 


Z2 


Wl + 





\ ^4 J 




\ 


) 




\ J 




I 1 / 



W2 



By the nature of this partial prolongation, the W2 vector must be of the form 
(0 1)^. The syst ems (x,u) and (z,w) are dynamically equivalent. Through 
a change of basis, transform the wi vector into (0 1)"^. 



/ 



-Z2 




1 



-zi 
1 









\ 


( ' \ 














Z2 







\ / 


V 1 / 
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This corresponds to the change of coordinates 

{zi,Z2,Z3,Z4) ^ { Z3 ~ Z1Z2, Z2, Z4, Zl). 

The change of coordinates is a static equivalence between (z.w) and (z,w), and so 
we have yet another system dynamically equivalent to (x, u). 



W2 



The (z,w) will be a partial prolongation of a three state system. In this case, 

Wi ^ Vi, W2 ^ V2- 

The numberings were chosen so that the final equations of the control system end 
up in this particularly nice form. 



/ 








-Z3Z4. \ 












Z2 






Z3 














Z3 









+ 





Wi + 


1 


V 


Z4 / 






) 




1 1 ) 




[0 J 





By construction, the systems (x, u) and (y, v) are dynamically equivalent. What 
the process does not tell us is if this equivalence is strictly dynamic, for it could 
easily be static as well. In this example, however, the (x, u) system is one of the 
classes of static equivalence given in Elkin [2\ and Wilkens [13j . while the (y,v) 
system is clearly static equivalent to a distinct class (y, v) 



(3.6) 

following the transformation 







( 




- 





Vz / 




^ 2/2 




vi = -y2Vi, 



1 = 1,2,3, 

V2 = V2. 



It is interesting to note that unlike the original system (3.5) in our equivalence 



(3.6 1 decouples into two smaller and separate systems: the first equation involves 
just yi, vi, while the other two equations involve only y2, ys, V2. This equivalence 
also converts a nonlinear system (x, u) into a linear one (y,v). Both decoupling of 
equations and linearity greatly simplify the analysis of solutions of control systems. 

Not only does the process presented above tell us that (x, u) and (y, v) are dynam- 
ically equivalent, but through some back tracking, it gives us the explicit equivalence 
maps. 



{Xi, X2, X3, Ui, U2, 



{yi, y2, VS, Vi, V2, ...)^{- 



Vi 



X1X2, U2, X2, 



2/3, 2/1 



-X1U2, U2, 

V1V2 - y2Vi 



■) 



2/2, 



■) 



ysvi 

y2 ' ' 2/2 ' 2/2" 

This simple example also shows why it is necessary to consider dynamic equiv- 
alence on open sets. In this case, we would need to restrict our equivalence to the 
open set where 1/2 7^ 0. 
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4. Previous Results 

The first theorem of this section is one of the most important, yet simplest to 
state, properties of dynamic equivalence. It can be found stated in a compatible 
form in [3_, but the following theorem and its proof, which are more in line with 
the terminology of this paper, can be found in [9_. 

Theorem 3. The number of control variables is an invariant of dynamic equiva- 
lence. 

Note that while this theorem states that dynamically equivalent systems must 
have the same number of control variables, they may have different numbers of 
state variables. This is most obviously illustrated by Theorem [2] A system with n 
states and s controls is equivalent to its prolongation, which has n + s states and s 
controls. Thus the number of states in a system dynamically equivalent to a given 
system is unbounded. 

Recall that a submanifold of an affine space is called ruled if, given any point 
of the submanifold, there is a line that passes through that point and that is con- 
tained completely within the submanifold. Classic examples of ruled submanifolds 
are planes, cylinders, and the hyperboloid of one sheet. We will abuse this termi- 
nology slightly and still call a submanifold ruled if it is the intersection of a ruled 
submanifold with a possibly bounded open set. A control system is called ruled if, 
when viewed as a subbundle E of the tangent bundle TAI, it defines at every point 
X a ruled submanifold Ej; of the tangent space T^M at that point. 

To state what is probably the most significant result to date in dynamic equiv- 
alence, some notation must be established. For j < k, let tt^j be the canonical 
projection from J*'{M) to J-'{M). Obviously TTfe.fc is the identity. For any open 
set n C J^{M), define the subset C j\M) by 

The following is due to Pomet [lOj. 

Theorem 4. (Pomet) Let E and A be control systems with state manifolds AI 
and N of dimension m and n, J, K two positive integers, and U C J'^^^{M), 
V C J^'^^{N) two open sets satisfying 

(4.1) Z^iHEc (Z^nE,7+i)i and Vi H A C (V n Ak+i)i. 

// E and A are dynamic equivalent over hi and V , then 

• if m > n, then E is ruled in lAi. 

• if n > m, then A is ruled in Vi. 

• if m = n, then 

— (real analytic case) iflAi H E and Vi n A are connected, either E and A 
are ruled inlAi and Vi, respectively, or they are locally static equivalent 
over lAi and Vi. 

— (C°° case) there are open sets R,SciUi and 7?,, 5 C Vi with 

(1) Ui = RU S = RU 

(2) Vi = 7^ U 5 = 7^ U 5, 

(3) E and A are ruled over R and TZ, 

(4) E and A are static equivalent over S and S. 
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The condition |4.1| basically says that nothing is lost when either prolonging the 
control system up or projecting the open set down in the jet spaces. In fact this 
containment is an equality; the reverse inclusion follows directly from the defini- 
tions. 

Recall that every system E is dynamically equivalent to its prolongation E2- 
Since the dimension of the state space of E2 is larger than the dimension of the 
state space of S, this theorem guarantees that E2 is ruled. Of course we already 
know that S2 is affine linear, so in this case the result is trivial. A natural question 
that arises from this, and one that is partially answered by this paper, is this: 

Given an affine linear control system, when is it the prolongation of a smaller 

system? 

At the moment, this question has not been answered in its full generality, here 
or elsewhere. In an attempt to partially address this issue, this paper will classify 
control systems of low dimension that are affine linear up to dynamic equivalence 
in Chapter [12] The methods used to do this rely on a previous classification of 
affine linear control systems under static equivalence. Gardner and Shadwick first 
classified control systems with two state variables and one control variable under 
static equivalence [5]. Wilkens then solved the problem for three states and two 
controls [12 . For the complete classification of affine linear systems under static 
equivalence with at most three states, which I present here without proof, see |2]. 

In the following theorem, n represents the number of state variables Xi, and Uj 
are control variables. Given a control system E : x = f(x,u) with state space M, 
we say that a point p € Af is regular if there is a neighborhood of p on which the 
rank of E, defined to be the rank of -i^, is constant. 



Theorem 5. (Elkin) An affine linear control system (3.4 1 with n < 3 states is 



locally static equivalent at a regular point p to one of the following systems: 

• n — 1 

Xl =0, Xl = 1, Xl — Til. 

• n = 2 



Xl 
X2 

Xl 

X2 



= 


f Xl 




Xl 




= ' 


[ X2 


: i 


±2 


= 


= Ul J 


Xl 


= Ul 


Xl 


= Ul 


= 1 ' 1 


X2 


= Xl ' ] 


±2 


= U2 



• 71 = 3 



Xl 


= 


r Xl 


^ 1 


Xl 


= Ul 


( Xl 


= Ul 


±2 


= , < 


1 


- , < 


X2 


, 


I X2 


= 1 , 


±3 


= 


[ ±3 


= 


> ^3 


= 


[ is 


= 


Xl 


= Ul 1 


Xl 


= Ul 


Xl 


= Ul 


Xl = 


Ul 


X2 


^ Xl A 


X2 


= Xl ,< 


X2 


= Xl ,< 


X2 = 


H{x)ui 


±3 


= [ 


±3 


= 1 


±3 


= X2 


. ^3 = 


l+a:2Ui 
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where H{x) is an arbitrary function with is nonzero 



Xi 


= Ul 


±1 = 


Ul 


Xi 


= Ml 




X2 


^ U2 , 




T • 1 


X2 


= U2 




±3 


= 






. ^3 


= M3 




ii 


= Ul 


±1 = 




±1 




Ml 


±2 


= M2 , < 


i2 ^ 


U2 , < 


X2 




M2 


±3 


= X2 


. ^3 = 


a;2Ui 


. ^3 


= U 


' X2U1 



The following theorem takes the problem of classifying control systems with one 
control variable under dynamic equivalence and reduces it to the simpler case of 
static equivalence. While this theorem has been known for some time (see [9] for 
one example) , a new proof of this theorem in the framework of this paper will be 
given in Chapter [8] 



Theorem 6. Let the control systems S, A in (3.1) he dynamically equivalent with 
5 = 1 control variable and m, n state variables, respectively. If m — n, then the 
systems are in fact static equivalent. If m < n (m > n), then the systems are static 
equivalent after a finite number of prolongations of the smaller system I] (K). 



5. The Equivalence Problem 

Given a manifold M, a framing on M is a collection {X^}"^]^ of smooth sections 
of the tangent bundle TM such that for every p g M, the collection of vectors 
{(Xi)p}"^j^, called a frame, forms a basis for TpM . A coframing is simply the 
dual of this notion, i.e. a collection of 1-forms {w"'}"^i (smooth sections of the 
cotangent bundle T*M) such that {(i^"')p}"=i forms a basis for T*M for every 
p € M. Every coframing u)^ has a corresponding framing Xi for which a;-'(Xi) = 51. 
A local framing/coframing is simply a framing/coframing defined on an open 
set U C M. 

An equivalence problem [4 can be stated in the following way: Let M" 
and N"" be smooth n-dimensional manifolds and G C Gi(n,M) a subgroup. Let 
LOu — {uJij, . . . , i^u)^ and ily — {^v, ■ ■ ■ , ^v)'^ t)e local coframings of U C AI and 
V G N, respectively, chosen in some geometrically natural way. We wish to find 
necessary and sufficient conditions that there exists a diffeomorphism (p : U ^ V 
such that 

where jvu ■ U ~> G. A common abuse of notation, one which will be used in this 
paper, is to drop the puUback from the notation where the map ip is clear from 
context: Qy = "fvu^u- 

For example, suppose we are given manifolds M and N with Riemannian metrics 
ds'^ and dS^, respectively. We can locally diagonalize the metrics on open sets 
U cM and V CN such that 

ds'^Y^i^hf. ds'^Y^^^vf- 

i i 

The problem then is to find necessary and sufficient conditions such that a diffeo- 
morphism ip : M ^ N exists such that ip*riY = ^yijuju, where ^yu is an element 
of the orthogonal group 0{n). 
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The goal of this paper is to adapt the framework of an equivalence problem 
to dynamic equivalence. Then, using methods of exterior differential systems, we 
will classify a collection of control systems. What makes the dynamic equivalence 
problem tricky is the unboundedness of the size of the potentially equivalent state 
manifold, and hence also the lack of diffeomorphisms. A diffeomorphism if : M ^ 
N cannot exist due to differences in dimension. In fact, strict dynamic equivalences 
are defined in terms of submersions rather than diffeomorphisms. This difficulty 
due to submersions persists through any finite number of prolongations. To solve 
this problem with submersions, in the next section we will simply make everything 
the same size: infinite. 



6. Infinite Prolongations 

The trick to dealing with our submersion woes is through prolongation, an idea 
introduced in section [3j Recall that a control system on M 

(6.1) S: = /,(x,u), l<i<n, 

can be represented by 

as a parametrization of S inside M x TM. A basis for the space X"*" is 

UJi — dxi — fi{x, u) dt, 1 < i < n. 

The forms uJi are the puUback to S by the inclusion map of the contact forms 
dxi — Xi dt on M. X TM, where M x TM has coordinates (t, Xi, Xi).Th.e collection of 
1-forms {dt,uj'' , duj} forms a coframing on E that encodes the information of the 
control system. 

Prolongation of ( |6.1| yields a system E2 given by the equations 

X = f(x, u) 
ii = u 

with state variables x, u and control variables u. Thus a suitable coframing on S2 
that encodes the information of the original control system and E2 is 

= dt, 

ujf = dxi — /i(x, u) dt, 1 < i < n, 

= duj ~ Uj dt, I < j < s, 
uj'j = duj, I < j < s. 

Define the infinite jet bundle J'°°{M) as the projective limit of the finite 
jet bundles J°°{M) = ^im J'^(A/), endowed with the projective limit topology. 

K 

Let Soo and Aoo be the projective limits of the prolongations of the control sys- 
tems E and A, respectively. By repeated iterations of the prolongation process 
above, a suitable choice for preferred coframings on Eco and Aoo with coordinates 
{t, X, u, li, ii, . . .) and (i, y, v, v, v, . . .), respectively, which encodes the information 
of the respective control systems is as follows. 
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(6.2) 



ft 



02 



dt 

dx — /(x, u)dt 
du — xidt 
du — iidt 



V 



dt 

dy- g{y,v)dt 
dv — wdt 
dv — wdt 



\ 
J 



The covectors a;* and Jl* are n-dimensional for z = and s dimensional for z > 0. 

Now we should take a closer look at what happens to the mappings involved 
in the definition of dynamic equivalence under this infinite prolongation process. 
Given a map $ : Sj+i — > N, as in the definition of dynamic equivalence, define the 
l^th prolongation of the map. denoted as the map that makes the following 
diagram commute on solutions. 



PJ+k+lX 




In other words, 



{pk o $) ( pj+ix{t) ) = ( pj+k+ixit) ) 

for solutions x{t) e M. 



Now define $o 



by $00 — '^'[fc] the obvious fashion, i.e. for 



projection the projection map tt^ that takes an infinite jet to the fc"^ jet, 

TTfc O = $[fc] O TTj+k + 1- 

Let ^E" : A/^ +1 — M be the map used in section |3] in the definition of dynamic 
equivalence, and define 4' 00 similarly. From the definitions of dynamic equivalence 
and prolongation, it is simple to show that 

* o ^[K+i] o pj+k+2 = Ida 
is the identity on curves in AI. Finite prolongation of this relation shows 

(6.3) *[fc] o $[A'+l + fc] o Pk,J+K+2+k = Idk 



is the identity on curves in S^. Taking the limit of (6.3 1 as k tends to infinity tells 
us that 



'I' o <t> —Id 
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is the identity on ^cxj- Similcirly 

is the identity on J'°°{N), and we can conclude that $00"^ = ^00 and that $00 is 
a diffeomorphism. 

To recap, in order to pose an equivalence problem for dynamic equivalence, we 
need a diffeomorphism between spaces. The problem with djmamic equivalence 
is that the maps used in the definition of the equivalence can never give us a 
diffeomorphism at any finite level (unless the equivalence is actually static). By 
passing to the infinite prolongation, the submersions become diffeomorphisms. We 
obtain the nice transformations we wanted, and now the issue is that we have to 
work on infinite-dimensional spaces. 



7. Group Action on the Infinite Prolongations 

Now that we have our diffeomorphism between infinite jet bundles, we would 
like to know the form of our group action G. Instead of working with a sub- 
group of GL{n, M), what wc have now is a group of transformations that take local 
coframings of T*T,ao to local coframings of T*Aao. In an equivalence problem of 
finite dimensional objects, (p*u} = 7$!, 7 is essentially the pointwise Jacobian of the 
diffeomorphisms tp. The same is true in the case of infinite prolongations. 

Suppose we have a transformation (t, x, u, u, . . .) i-^ (i, y, v, v, . . .) such that 
t^t. Suppose y = y(x, u, li, . . . , u^"''), i.e. yu(^) is nonzero and yu(fc) = for all 
k > J. If V = v(x, u, li, . . . , u^''^^), we need to know first of all how Ji is related to 
J. 

On the one hand, we can directly compute the time derivative of y using the 
chain rule. 

^y(x, u,u, . . . ,u("')) 

yx(x, u, u, . . . , u('^))f (x, u) + yu(x, u, li, . . . , u'^'^^)u 
+ ... y„u,(x,u,u,...,u(^))u^^+'^ 

On the other hand, y = g(y, v). 

^ = g(y(x,u,u,...,u(^)),v(x,u,u,...,u(-^i))) 

Comparing these two versions of ^ shows that v = v(x, u, li, . . . , u('^+^'). Thus 
we have the following theorem. 

Theorem 7. yu(j) is nonzero and y^m = for all k > J if and only if v^(j+i) is 
nonzero and v^^(k) = for all k > J + 1. 

This relation and its repeated derivatives with respect to t show that 




v« =vW(x,u,u,...,u(-^+'+i'). 
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TheoremjTjrelates to our coframing as follows. Here we are omitting the puUbacks 
from our notation. 



dy - g(y,v)dt 



dy = d(y(x,u,u, ...,u(-^))) 

i— 



i=0 



X, u, . . . , u^-^+^M dt 



AS](dx - f (x, M)dt) + ^ A°+i (dvS^ - u'-'+^^dt 



i=0 



where A°, < j < J + 1, are matrices of functions of x, u 



,u('^+i). The fact that 
dynamic equivalence is time independent and takes solutions to solutions implies 



that there is no additional A^^j^ dt here. 



Similar calculations for dv^*) — v*^'+^^(it imply that our preferred coframings (6.2) 
transform in the following way 

. \ -1 



(7.1) ^lo^ = Ai, 

where Jl,A,ci; have the form 



V ■■ ) 



,J+1 

,J+2 



1 


Olxn 


Olxs 


Olxs 


Olxs 


Olxs 


Olxs 


Olxs 


Onxl 


A" 


A? 




AO 


Osx s 


Osx s 


Osxs 


Osxl 


Al 


A} 


Al 


■ Al 


Al 


Osx s 


Osxs 


Osxl 


A2 


A? 


A^ 

-^2 


• A2 


A^ 


A^ 


Osxs 



V 

and the A* are submatrices of the following sizes. 



matrix 


AO AO Aj, A} (z,j>l) 


size 


n X n n X s s x n s x s 



A matrix A of the above form for a fixed J may have an inverse matrix similar 
to the above form with arbitrarily large K. For example, composition of dynamic 
equivalence maps leads to arbitrarily large J and K. 

From here on out, for any statement or theorem about A, an analogous statement 
or theorem also holds for A~i unless otherwise noted. These have been omitted 
for brevity. Submatrices of A (A^i) will be denoted by uppercase A* (lowercase 
a*). Individual entries of these submatrices will denoted by {Aj)f ((a*)f), which 
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are functions and thus not bolded. If a particular submatrix is in fact a scalar, 
which happens when s = 1, then no bold face type will be used: A^ . 

Theorem 8. Given a dynamic equivalence $J^r2 = Aa; with adapted coframings 



= V alli>\. 



Proof. This proof is by induction on i. The case of i = 1 is obvious. For « > 1, 
consider d(r2*). Where an equivalence sign = is present below, it is because we are 
considering the equation modulo the linear span of {w", . . . , a;'^+'+^}. Keep in mind 
that we are working with vector equations here. Recall that a;" = dx — f(x, u)(it 
is n X 1, and ui^ = du^^-i) - vS^Ut, Q.' = dv(-'-i) - v^^^dt are s x 1 for j > 1. 
It is straightforward to verify in coordinates that du!^ = — a;^+^ A dt and dW = 
-rjJ+i A dt for j > 1. 



On the one hand, 



= -dv(*) A dt 
= -Jl'+i A dt 



-A^^,u;'^+»+2 A dt. 



On the other hand. 



J+i+l 



^ [d(A})Aa;^+A;.d(a;^)] 

J+i+l 

^ [d(A}) A a;^' - A}a;^+i A dt] 



^ -AV.+i^''+^+'Adi. 

Since the a;-' form a coframing, they are linearly independent. Thus we can conclude 
that 



□ 



While this does not completely characterize the group action of dynamic equiva- 
lence, it will be sufficient to prove a result in the next section that classifies dynamic 
equivalence in the case of one control variable. Later sections will narrow down what 
this group A looks like; however, we will never completely characterize it. What 
we do prove about A will be sufficient for some non-existence results. 

8. Scalar Control 

The following theorem about dynamic equivalence in the case of one control 
variable has been known for some time. What is presented here is a proof based on 
Pomet's work [9] that has been adapted to this framework of coframings on infinite 
jet bundles. It reduces all dynamic equivalences of control systems with just one 
control variable to the case of static equivalence. 
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Theorem 9. Let the control systems 'E, A in (3.1) be dynamically equivalent with 
5 = 1 control variable and m, n state variables, respectively. If m — n, then the 
systems are in fact static equivalent. Ifm < n, then the systems are static equivalent 
after a finite number of prolongations of the smaller system E. 

Proof. Let A = (Ap and A^^ = (ap as before. 

If TO < rt, prolong S until m — n. Suppose the coframings of E, A in (3.1 1 pull 
back as in (7.1). Suppose there exist nonnegative integers J and K such that x^(j) 



and Yy^iK) are nonzero. In Theorem 11 in the next section, it is shown that it is 
not possible for just one of J or iiT to be —1, i.e. A* = for all j > i if and only 
if a* = for all j > i. So both J and K must be nonnegative for a strict dynamic 
equivalence to exist. 

By the computations in the previous section, Aj^-^ is a nonzero n x 1 matrix. 
Likewise, a)^^^_^^ is a nonzero function for all i > 1. Because Aj^j^ is a nonzero 
n X 1 vector, and a^^^j_|_2 is a nonzero function, their product A^j^-^^a'^^j^^ is a 
nonzero n x 1 vector. However AA~^ is the identity. Therefore Aj^j^a^"^^j^2; 
which is an off diagonal n x 1 entry since < K + J + 2, must be an all zero n x 1 
vector. 

This is a contradiction. Thus J and K cannot exist, and x^(j) = y^(K) = for 
all J,K>0. This shows that the equivalence is in fact static. □ 



9. Group Adaptations for Two Controls 

The last section dealt with the case of a scalar control, in which dynamic and 
static equivalence are one and the same. Now we will work on the next simplest case 
of two controls (s = 2) with J = K = 0. In the case of one control variable, there is 
essentially no "room for freedom" to allow a true dynamic equivalence, aside from 
prolongations. With two control variables, there is now "room" to have a strict 
dynamic equivalence, but just barely. While larger values of J and K increase the 
flexibility of possible dynamic equivalences, in this section we will show that there 
is really only one way to have a strict dynamic equivalence of two systems with 
J = K = 0. 



9.1. Nonautonomous Static Equivalence. Recall the notation we have devel- 
oped thus far for the puUbacks of our preferred coframings. Note the equivalent 
submatrices A.l_^_^, i> 1, from Theorem 8 



V : J 



1 


Olxn 


0lx2 


Olx2 


0lx2 


Oi> 


<2 


Onxl 


A" 


A? 


Onx2 


0nx2 


o„, 


<2 


02x1 


Ai 


A} 


Ai 


02x2 


02> 


<2 


02x1 


A2 


A? 


A2 

^2 


Ai 

-^2 


02> 


<2 



, ,0 
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1 




Olxn 


0lx2 


0lx2 


Olx2 


Oi> 


<2 




1 




a? 


0nx2 


Onx2 


o„> 


<2 




1 




a} 


4 


02x2 


02> 


<2 


02X 


1 


4 


a? 


ai 


4 


02> 


<2 



V 



V ■■ J 



In what follows, we will refer to a group element G, 



G = 



V 



1 

Onxl 
02x1 
02x1 



Olxn 
g§ 



0lx2 
0nx2 

g! 
g? 



Olx2 
Onx2 
02x2 

gi 



0lx2 
0nx2 
02x2 
02x2 



Olx2 
Onx2 
02x2 
02x2 



that acts on our coframings as nonautonomous static equivalence, meaning 
G* = for all i < j. This terminology arises from the fact that such G arise as 
the Jacobian of a time-dependent static equivalence x = x(x, t) on the contact 
system of the infinite prolongation Eoo- Unlike the matrix representing a true 
static equivalence, G allows changes of variables such as Xi i-^ Xi + t. Note that 
such equivalences take a coframing on Eqo to another coframing on Sqo (or Aqo to 

Aoo). 

As in the case of dynamic equivalence, we wish to require that the following 
structure equations are preserved by nonautonomous static equivalence. 

(9.1) 



dn' e span A n^'} mod n\ < j < i. 



This additional condition allows us to simplify the form of G much like we did for 
A in Theorem [8j The proof is identical to that of Theorem [8] with J = — 1. 

Theorem 10. Given a nonautonomous static equivalence G that preserves the 
structure equations (9.1 1, g* — g\ for all i>l. 



A straightforward calculation in coordinates shows that every static equivalence 
is a nonautonomous static equivalence, but of course the converse is not true. 

Later we will be showing that every dynamic equivalence with J = K — can be 
factored into a constant matrix composed with nonautonomous static equivalences. 
This result will be key in proving the main classification results of this paper. 

Theorem 11. A is a nonautonomous static equivalence, i.e. A* = for all i < j, 

if and only if A^^ is also a nonautonomous static equivalence. 

Proof. If Ai = 0, then Aq is a rank n matrix, hence invertible. Let the submatrices 
of A^^ be denoted by a*. Since AA^^ = Id, the off diagonal element Aq&I must 
be zero. Because Aq is invertible, this means a° = 0. By Theorem^ a° is zero if 
and only if is too. Theorem [s] completes the proof since al^i — a| for all z > 1 
and a* = for all j > i. Therefore A~^ is a nonautonomous static equivalence. □ 

9.2. Factoring A. In the following section, we will prove several theorems about 
the rank of certain submatrices of A. This chapter will culminate in the final 



theorem, theorem (14 1, which states that we can factor our dynamic equivalence in 



a special way: A = gSQ. The G and G are two nonautonomous static equivalences 
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which encapsulate the traditional change of variables, as in static equivalence. The 
iS is a fixed constant orthogonal matrix which incorporates the mixing of higher 
derivatives into dynamic equivalence. 



Theorem 12. Given a strictly dynamic equivalence A with s 
A2 (a 2 X 2 submatrix) has rank 1. 



2 and J = K = 0, 



11 



Assume the rank 



Proof. We know that Aj cannot have rank zero by Theorem 
of A2 is two. Then through a change of coframing a> = Guj via static equivalence 
G, it can be arranged that the elements of A = AG^^ look as follows. 



V : J 



( 



V 



1 

Onxl 
O2XI 
O2XI 



Olxn 


0lx2 


Olx 


2 


Olx 


2 


Oi> 


<2 


A° 


A? 


On> 


2 


On> 


2 


o„> 


<2 


02x2 


02x2 




X2 


O2X 


2 


02> 


<2 


02x2 


02x2 


02X 


2 


/rf2 


X2 


02> 


<2 



\ ■ J 



We have ^^^2^ = a>^^^ for j < 1. By the nature of puUbacks, this also means 
) u^^^ — W . However this means that A^^ now looks as follows. 



^2 



1 


Olxn 


Olx2 


0lx2 


Oi> 


<2 


0i> 


<2 


nx 1 




a? 


0nx2 


0„) 


<2 


0„> 


<2 


2x1 




a} 




02> 


<2 


02> 


<2 


2x1 


02x2 


Id2x2 


02x2 


02> 


<2 


02> 


<2 



V : / V : J 

In particular, = a| = aj. By the above argument this means that 
equivalence is static. This contradicts J — K = 0. Therefore the rank of Aj must 
be one. □ 



V J 

and the 



Theorem 13. Given a dynamic equivalence A with s = 2 and J — K — 0, Aj (an 

n X 2 submatrix) has rank 1. 

Proof. The rank of Aj is either 0, 1, or 2. If the rank is zero, then the equivalence 
is static. Consider (AA~^)2 = A^aj = 0„x2- If rank of A\ is two, then A° has a 
2 X n left inverse, and we conclude a.2 — 02x2- However this again implies a static 
equivalence, so the rank is not two. □ 



The plan now is to use this knowledge of the ranks to normalize A via non- 
autonomous static group actions to uj and fl. This will isolate the dynamic part 
of the mapping to one very specific form Cl = Slj, where Cl = G~^f2, u) = Gw, lj 
and Jl are our preferred coframings (to be determined later), and G, G are non- 
autonomous static group elements. An explicit example of how this is done will 
follow in the next section. 

Starting with the fact that A'j' has rank one, we know it can be normalized 
to the following form through Gauss-Jordan elimination, which in this context is 



20 



MATTHEW W. STACKPOLE 



nonautonomous static equivalences applied to the coframings w and fl. 

/ 1 \ 



A? = 







Recall that all = {ujj) and fl^ = (fij) are vectors. If we add multiples of ui^ to 
col, ^6 can eliminate the first row of Aq. Note that this can be accomplished by a 
static group action. 

/O ••• 0\ 

Aq — 

y * • • • * y 

Since the n x (n + 2) matrix (Aq | A°) must have rank n for A to be invertible, the 
last n — 1 rows of Aj] must have rank n — 1. This allows us to normalize the rest 
of Aq via a static group action. 



AS 

^0 



Olx(n-l) 



*(n-l)x(n-l) 

The first n + 1 rows of A have now been reduced to ones and zeros. 

Since the rank of A^_|_j^ is one, there exists a non-autonomous static equivalences 
applied to both coframings w and Ct that yields a new coframing with 

A^ -("O 
A^+i - I 1 



can be eliminated by a nonau- 
In fact anything to the left of 



Everything to the left of the ones in each A]^-^ 
tonomous static equivalence that redefines ujIj^i- 
or below a one in the matrix A can essentially be absorbed by a non-autoiioinous 
static equivalence that redefines either uj (horizontal zeros) or S7 (vertical zeros). 
For ease of notation, these newly redefined coframings, which differ from the orig- 
inal preferred coframings by non- autonomous static equivalences, will be still be 
denoted with and u). This leaves the following simplified form of A. 

/lO 00000000 ■•■\ 



V 



1 





























■•■ 














1 

















■■■ 








Id(n-l)x(n-l) 























■■• 





(Ai)i 





(Ai)i 




















••• 




















1 











■■■ 





(A^)i 





(A^)i 





(A^)i 














••• 


























1 





••• 





(A^)i 





(Af)i 





(A^)l 





iM)i 








••• 
































••• 



/ 



Now if {Aq)\ is zero, one of the other (Aq)}, i > 1, must be nonzero. This follows 
from the fact that A^-^A = Id, in particular ((A~^A)q)J = 1. If (Aq)| is zero, 
there is an i > 1 such that (a°)}(AQ)J is nonzero. But (a'')i[ being nonzero implies 
K + 1 > i > 1. Since we are restricting our consideration to K = 0, this cannot 
happen. Therefore (Aj)} must be nonzero. Since (Aj)} is nonzero, it can be scaled 



DYNAMIC EQUIVALENCE OF CONTROL SYSTEMS VIA INFINITE PROLONGATION 21 

to unity through a nonautonomous static group action. All of the other (Ag)} can 
then be eliminated through non-autonomous static group actions (adding multiples 
of rows in this case) . 

It can similarly be shown that when J ~ K ^ 0, (A^^^)} is nonzero and can be 
scaled to unity. All entries below them can be made zero. By examining A^^ A = Id 
one can also check that any of the (A^)} being nonzero leads to + 1 > 2, and 
therefore (A*)} — 0. 

Finally all the group freedom of A has been utilized through non-autonomous 
static group actions on uj and O, and what is left is the following constant matrix. 



(9.2) S 



1 


Olxn 













Onxl 


Oix{n-l) 
-f'i(n-l)x(n-l) 


1 

0{n-l)xl 0(„_i)xl 


Onxl Ojixl 


Onxl Onxl 








1 ^lx(n-l) 
Oix{n-l) 







1 












Olxn 
Olxn 


1 








1 

















V 



J 



It is easy to check that S is orthogonal, i.e. S ^ — . Thus we have proved 
the following theorem. 



Theorem 14. Given preferred coframings and ft (6.2 1 on Eqo md Aqo 

control systems E and A, respectively, with s — 2 and a dynamic equivalence <i>, 
with J = K = taking Eqo to Aqo, the coframing pulls back as follows: 



for 



G S G u) 



where G and G are nonautonomous static equivalences and S is given by (9.2) 
above. 

This theorem means that, up to nonautonomous static equivalence, a dynamic 
equivalence with J = K = has a very specific form which is encoded in this specific 
orthogonal matrix S. Most of the apparent complexity of dynamic equivalence 
actually arises from static equivalence on either side, and the essence of dynamic 
equivalence is actually quite simple. 



10. FACTORING THE DYNAMIC EQUIVALENCE: An EXAMPLE 

What we have shown so far is that, given a dynamic equivalence = Au; where 
J — K — and s = 2 (n is still arbitrary), we can decompose the group action 



A = gSG where S is defined by (9.2 1 and G and G are non-autonomous static 
equivalent group elements, i.e. 



/ 
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Recall that nonautonomous static equivalence is not a true static equivalence 

(x,u) (y(x),v(x,u)). 

Unlike static equivalence, which is autonomous (time- independent), nonautonomous 
static equivalence can have explicit time dependence, for example, i— > a;, + t. 
Its group action does not preserve the algebraic ideal {dx}, just the ideal {dx — 
f(x,u) dt}. This equivalence is more general than static equivalence. 

Let us phrase the problem now as follows. Dynamic equivalence looks like = 
Au where A = gSG. We can attack this problem in steps. First we will consider 
the coframing ft = SGui. Then what remains will be the non-autonomous static 
problem ft = gCl. 

ExEimple 4. 

Let us consider the following two dynamically equivalent systems: 

xi =ui yi= vi 

X2 = U2 2/2 = V2 

±3 = X2U1 1/3 = y2 

An actual dynamic equivalence is given by the following maps between the infinite 
jet bundles. 

$co(x, U, U, ...) = {X1X2-X3, U2, X2, X1U2, U2, ■■■) 

^^^(y>v,v,...) = ( Ui/t/2, ys, 2/3^'l/y2 - 2/1, ■••) 

Here is a coframing for each of the infinite jet bundles. The choice of ^3, while not 

obvious, is not arbitrary. We will see why in a later section. For this example only 
the dt piece of the coframing has been left out. Since t ^ t, this would just add a 
one and many zeros to the matrices. 





dxi 


— Uidt 




dx2 


— U2dt 


dx3 


— X2Uidt 


- X2{dXi 


dui 


— uidt \ 




du2 


- U2d.t J 




dill 


— uidt \ 




du2 


— U2dt ) 





dyi — vidt 
dy2 - V2dt 
dys - 2/2 rfi 
dvi — vidt 
dv2 — V2dt 
dvi — Vidt 
dV2 — V2dt 



The pullback of $00 is straightforward to calculate. For the rest of this section 
the pullback notation will be suppressed in order to emphasize and clarify the 
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methods being used. 

= dyi — vidt 

= d{xiX2 - X3) - (X1U2) dt 

= X2 dx\ + x\ dx2 — dx3 — {X1U2) dt — {x2Ui) dt + {x2Ui) dt 
= — [dx-i — x2U\dt — X2{dxi — uidt)] + x\ {dx2 — U2 dt) 
= -W3 + X1CO2 

O2 = dy2 — V2dt 

= dU2 — U2 dt 
= ^2 

^3 = - y2dt 

= dX2 — U2dt 
= 0J°2 

Q\ = dvi — vi dt 

= d{xiU2) - {uiU2 + X1U2) dt 

= U2 {dxi — ui dt) + x\ (du2 — U2 dt) 
= U2CJ1 + X-1UJ2 

fll = dv2 — V2 dt 

= dU2 — U2 dt 

= w| 

The puUback put in matrix form looks as follows. 







f 





Xl 


-1 

















• 


•• \ 


^2 




















1 











• 




eg 


0° 









1 




















• 




col 


vl{ 






U2 











Xl 











• 




ujI 


n\ 
























1 





• 




coi 


vl{ 






U2 








U2 


Ui 





Xl 





• 




ujI 


nl 






























1 • 




, ,a 
^1 


\ ■■ J 




[ 




















/ 


W2 



We will now follow the algorithm for producing S. This amounts to a series of 

row or column operations which arc static equivalences on the Q, or uj coframcs 
respectively. We will use the notation of a typical introduction to linear algebra 
course to represent these operations, i.e. i?2 R2 + R3 means to replace row 2 with 
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row 2 plus row 3. Note that not every row operation is a legal static equivalence. 
For example, i?i i?i + R4 amounts to a; i-^- a; + u, which is dynamic, not static. 
First perform the following operations: 



i?l — y i?3 — > R2 — ^ Rl 

which results in the following coframing. 



/ ^"2 \ 




f 














1 











• 


•• \ 















1 




















• 






eg 











Xi 


-1 

















• 








n{ 






U2 











Xi 











• 






^2 


fii 
























1 





■ 






ui 


ii'i 






U2 








U2 







Xl 





• 






w| 
































1 • 






/ 1^* 


V ■■ ) 
























J 




/ ,3 

W2 



Next perform the operation 



to get this new coframing. 



R3 X1R2 - R3 















1 











••• ^ 





1 




















••• 








1 

















••• 


U2 











Xl 











••• 




















1 





••• 


U2 








U2 


Ul 





Xl 





••• 


























1 ••• 

J 



The first three rows of the transformation now look like the first three 
Continue by letting 

i?4 — XiRi 



rows of <S. 



i?4 



U2 



to yield the coframing below. 



xiVlI - nl 
{n\ - xinii) /u2 















1 











0... 





1 




















... 








1 

















... 


1 























... 




















1 





... 


U2 








U2 


Ul 





Xl 





... 


























1 • • • 



, .0 



V : / 
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Re 



Now the first five rows match S. One more operation 

Re - (U2-R4 + uiRi + X1R5) 

puts the coframing in the following form 

xif7§ - nl 



{nl - xinii) /u2 
^ 

nl - {U2 {U20.\ + Xi^l) + UiVtl + XxVi^ j /U2 



J 



/ 
1 
1 


1 









••• 

••• 
••• 


1 









1 


••• 

••• 





V : : : 


1 







••• 
1 ••• 



^,0 



V : ) 

and all visible rows now match those of <S. Continuing this process ad infinitum 
gives us new coframings that transform via <S. At present, this transformation looks 
like 

where G is the identity Id. To put it in the desired form, we simply invert the 
action on the left hand side. This results in the following factored transformation. 



n 






Xi 


-1 














1 

















■■■ 





1 














••• 


Xi 
















••• 














1 





••• 


—X1U2 — Ui 








U2U2 




U2 


••• 




















1 ••• 



S Id u) 



Note that this decomposition using non-autonomous group elements is not unique, 
however it was chosen so that the second non-autonomous group element of the 
latter equation was particularly simple (the identity in this case). Any problem 
with three states and two controls can be simplified in a similar way, as we will see 
below. 
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11. Three States and Two Controls 



11.1. Preferred Structure Equations. In the method of equivalence, described 
in Chapter [5] one important step is to work with an initial, preferred coframing that 
encapsulates the problem at hand and satisfies some particularly nice relations that 
ought to be preserved by the equivalence in question. In this section we will make 



one final refinement to our coframings (6.2) so that they satisfy some particularly 




nice structure equations that ought to be preserved by dynamic equivalence. 

Note that for a control system x = f(x, u) with n state variables and s < n 
control variables, the vector 



f(x,u) = 



must have rank^ = s. Therefore, by the implicit function theorem, a static 
equivalence always exists so that the above system is equivalent to x = f(x, u) 
where 

/ u^ \ 

Us 

/s+i(x,u) 

where xi — Xi for 1 < i < n up to reordering and Uj ~ fji^, u) for 1 < j < s. 

We will now, and for the rest of the paper, concern ourselves with the case of three 
state variables and two control variables. The above adaptation suggests altering 



f(i,u) 



(6.2 1 for the case of three states and two controls to the following coframing. 



-1 


= dt 






= dxi 


— ui dt 




dx2 


— U2 dt 




= dxs 


- /(x,u) dt 




= dui 


— ill dt 


uil 


— du2 


— U2 dt 



Here /(x, u) is a scalar function. Note that in this coframing, dw* = — o)*^^ Aw^^ 
for i > and j = 1,2. The outlier in this nice pattern of exterior derivatives is, of 
course, 

2 



3 

E. 



1=1 



With one more adaptation of the coframing, we can make even this structure equa- 
tion easier to work with. Let the following be our preferred coframing for the case 
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of n = 3 state variables, s = 2 control variables. 





= dt 




, lO 


— dxi 


— ui dt 


/ lO 

W2 


= dX2 


— U2 dt 


col 


= dx3 


-fdt^ 


col 


= dui 


— ill dt 


col 


= du2 


— U2 dt 



(11.1) = dx^- f dt~ fu^{dxi- ui dt) - fu2{dx2- U2 dt) 



Note that this cofraniing satisfies some particularly nice structure equations. 

duj^l — —uo\ A uo~^ 

doj2 — —^2 ^ 

doj^ = mod a;0 

dcol = -Lol'^^ A {j > 0, fc = 1, 2) 

We will take this coframing, along with the analogous coframing Jl in (y, v) coor- 
dinates, as our starting point. Let O — Glj and Cl = SO so that ft = Gft. In 
addition, we will require that at every step of our transformation of the coframes, 
LJ, tt preserves the following nice properties of the structure equations and their 
algebraic ideals: 



dco^ = — A LO 



-1 



mod Lij 







(11.2) 

dool = -Lol'^^ A cj-i mod {uj' \ 0<i<j }, {j > 0, k = 1, 2). 



dco^ = 



11.2. Reducing G. Consider the coframing ft = SGu). Since we plan on applying 
a generic G in the non-autonomous problem = gCI, G does not have to be 
completely generic. It can be simplified to remove some redundancies. For example, 
since cjg 1— >■ jig under S, there is no need to add an arbitrary multiple of to any 
other form through G since this can be taken care of with G. What follows will 
illustrate this more explicitly. 

We have coframings u) = Guj and ft = Sco = SGu). Recall that — G\ for all 
^ > 1 by Theorem [8] Consider the following identities. 

f]? = {Gl,)jcoUiGl)lu:'> + iGl,)lu;l + iGl)lcol + iGl)lu;l 

m = (GO)? + (GS)^ -2° + (GS)^ -3° 



n° = {G'o)lcoU{G"o)luj'2 + {G",)l 



nl = (gO)1c.o + (gS!)^^2° + (gS?)J^3° 

fll = (G^)? .0 + (G2)2 ^0 + (G^)2 ^0 + (^2)2 ^1 ^ (G?)i U^l 



{GDI + (Gi: 



2 



nl = {GDI coU (GDI 00^2 + {GDI 00^ + {GDI col + (GDI col 

nl = (Gil)? .0 + (G3)2 ^0 + ^0 + (G3)2 ,1 + (G3)2 ^1 

+ {Gl)l col + {GDI -2^ + (G})? co! + {GDl col 
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+ {Gl)l u^l + iGDlul 

+ {Gt)l Ji + {Gt)\ + {Gt)\ cjf + {GDI ujI + [GDI + {GDI u\ 



Now G will add arbitrary multiples of 0° and fig to every other part of the coframing 

in order to get the final coframing fJ. Since they are linearly independent, they do 
not need to be completely arbitrary We will not lose anything by letting (Gq)| = 
{G%)\ ~ 1 and {G^\ = {Gq)^ = 0. In fact all of the other terms above involving 
^2 and cjg may as well be set to zero since G will take care of these through 
nonautonomous static equivalence. 



= {GDI a;? + {GDI '^i + {GDI + {GDI + {GDI 
Cll = {Gl)lu;U{GDlu;l + {GDlu;'2 

Oi = {GDI a;? + {GDI u;l + {GDI c.^ + (G^)? + {G^l a;f 
+ {GDi ojI + {GDI 4 

= {GD\ujU{GDl^l + {GDl^l + {GD>! + {GDli^l 
m = {Gt>)i^U{GDl-^l + {Gt)lu.l + {Gt)lu:! + {Gt)lu:l 
+ {Gt)i + {Gt)l u:l + {GDI ^\ + {GDI ^\ 



Of course we are keeping careful note that every group reduction we have made is 
allowed due to the freedom we have in choosing G. 

Now it is clear that we may as well choose (Gq)} = 1, and thus we may also set 
any term involving lJ{ below f2} to zero since G will be adding arbitrary multiples 



{GDI 
{GDI 
{GDI 




of fi} to these. 



9l\ = {GDlu;U{GDiu^l + {GDliol 
fio = {GDlu;U^'2 
m = {GDl^'i+^'s 




{GDI u;l + {GDI coh + {GDI ^? + {GDI 



= {GDI ^1 + {GDI 
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nl = (Gt)?..i + (Gf)2a.2^ + (G|)f a;? + (G|)ia;| + (Gt)?a;? + (G^)ia;2^ 
+ {Gl)l Lot + iGDl Lol 

One entry in every tfj can be scaled to unity. Note that G\ is an invertible 2x2 
matrix, so that either the pair (G})}, (G})! or {Gl)l, {G\)l is nonzero. If the 
former pair is zero, then G would allow us to switch the roles of every Ci[ and 
for i > 1. Thus without loss of generality we can let {G\)l = (G})! = 1. The 
arbitrariness of G will then let us cancel out any terms below these scaled terms. 
For example, adding multiples of fi? and d\ to will get rid of the ui^ term in 
all the fj', i > 1. We can also scale the a;| term in to unity, and thus every u;| 
below can be eliminated. After this process of scaling one term per ^Ij and using 
this to eliminate the appropriate terms below, we are left with the following. 

fll = (Gj)f a;? + (Gl)f a;i+a;| 

rso _ (r'^y-)^ , ,0 I , ,0 

Clt = {G°)lco°+uj° 

nl = 

nl = {Gl)lojl + iG\)lu;l+ojl 

= ujI 

ni = {Gl)lu;! + {G\)lwl+u;l 

Cll = ojI 

nl = {Gt)lujl + {Gi)lu;t+ujl 



After all such redundancies are removed, this is what our group element, now 
called Q, looks like. 



1 
































••• ^ 





1 





























••• 







1 


























••• 





(G8)? 





1 























••• 














1 




















••• 





(GDI 








(GDI 


1 

















••• 


































••• 














(GDI 





(g\)1 


1 











••• 


































••• 




















{GDI 





{g\)1 


1 





••• 
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We will need to use the fact that a> is a coframing. Therefore exterior derivatives 
of the entries of Q can be written as linear combinations of these. Note that as far 
as we know, every d{Gj)f could be linear combinations of u)"^ for some unknown r. 
We will employ the following notation: 



(G})ti--^+EE(G})t>: 



We will show below that r is not arbitrarily large by looking at structure equations. 

By investigating dfl, we can further reduce the entries of Q. Until stated oth- 
erwise, the following equivalences = are modulo Q?, i = 1,2, 3. We will start with 



dill 



duj° + d{G°)l A w° + (G°)? duj° 

= (^U2fuiX2 ~ fxi — (Gq)^ + Ulfuixi + fulXsf + /a;3(Go)^ — fxsfui 

+ UlfuiUi + U2fuiU2^ ~ {GQ)l(^U2fu2X2 + Ulfu2Xi + U2fu2U2 ~ fx2 

+ /U2X3/ + Wl/«1«2 ~ /x3/«2) f^l A fj ^+ (^/ttiM2 — /«2«2 ('^1)1) ('^0)1 

~ ((<j'o)ii2(^i)i ~ fuiU2{Gl)l + fului — {Gq)\'^^ of A Q\ 
- {GDI Cll A 

Anything above that is not a multiple of Cll A or ^2 A must have zero 
coefficient. Of greatest interest at the moment is the term Clf A Since this 

cannot be here, its coefficient must be zero. 



(11.3) 



(G^)t = 



O^ 



There is also a Of A 0} term which must vanish. Through the above equation, this 
simplifies to the following. 

(^fuiU2 ~ /«2M2(^l)l) (^o)l + (/mi«2(^i)i ~ = 

Moving on, we will look at dfl2- 

= a;° + (G°)fc.? 
= dio^ + d(Gg)f A + (Gg)f dw? 

= {{GDI - (GDU) nl A + {{GDI - (GDI) A fl-' 

oo 

+ E ((Go)?;I - (gXUgDI - (G°)f;t,(q+i)f) 01+1 A 01 

i=l 

oo 

+ E(Go)?;'+i^^2 A 01 

i=l 

Similarly here it is the vanishing of the Of A O"! term that tells us 

(GDI = {GDI 
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The vanishing of the terms in the final two summations tells us 



(<^o)li2(^o)l> 

0, 




for alH > 2. We knew that 



had to be a finite sum, and now we have a bound on where that sum must terminate 

(«=!)• 

Now consider dfij. 



{GDI + u;l+u;l 

d{Gl)l A + {GDI rfw? + d{G°o)l A ojI + (C?g)f dul + (kvl 
\{Gl)li - {Gl)li){Gl)l + {GDli - {GDll - {Gl)li{Gl)l 



{G'ofdiGDi 



(Gj)2_^f7}AQ-^-O^AQ- 



{GDlpl A + [{GDI - {GDI - {GDI-,) ^1 ^ 

oo 

4 = 2 



T.(^DI'-+Ma^I 

4=2 



The relations that come from this calculation are these for i > 2. 



(GDI 

(gDI'.i 
(GDI- 

(GDia 



{GDi + (^"0)1,-1 

(GDliiGDl + iGDlD 



/^0\2,2/^0\2 /^0\2,2/^1n2 



Therefore we have found a bound on the sum for d(GD^'^ as well. 

Continuing this process for higher order terms yields the following important 
result 



{GiX\)l = {GDl + {i + mDl-i 



for i>l. 
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To review, Q now has the following form. 

(11.4) g = 



I 



1 


00 




















00 ••• 





1 00 




















00 ••• 





(GS])?10 




















00 ••• 





01 




















00 ••• 





00 


1 

















00 ••• 







(Gg)? 


1 














00 ••• 





00 








1 











00 ••• 





00 




-lO 


{G%)\ 


1 








00 ••• 





00 














1 





00 ••• 





00 










-lO 




00 ••• 



What we have boiled the problem down to now is the equivalence 17 = Gil, 
where the coframing contains three functions /, {G%)\, and (G}^\. 

Remark: An important but subtle point to take note of is the following: we have 
singled out through S as the piece of the coframing A° that contains higher order 
terms in w, and we have also singled out ^ by choosing an adapted coframing with 
df2§ = mod $7", and these two choices are compatible. 

This fact is actually quite easy to see. In our coframings, note that fjg = 0^3. 
Since G preserves the span of { fi?, fijj ^3 }) '^^s must be in the span of { ^^2) ^3 }• 
Thus cjg, which has the property that duJ^ = mod w". does not also get bumped 
up in the dynamically equivalent coframing to a higher order term. 



12. Dynamic Equivalence of Control Affine Systems 

Keep in mind at this point that we are concerned with dynamic equivalence, 
which is a weaker equivalence than static equivalence. The static equivalence case 
was dealt with first in the control linear case of three states and two controls 
by Wilkens and later by Elkin in the affine linear case up to four states. The 
representatives of the five distinct static equivalent control affine systems with three 
states and two controls put forth by Elkin are these: 

X\ = U\ 

(12.1) ±2 = U2 

±3 = /(x, U) 

where /(x, u) is one of the five following functions: 


1 

< X2 
X2U1 
I+X2U1 

In this section, we will finally put to use our previous results involving infinite 

prolongations and the factorization of coframing puUbacks. We show, using argu- 
ments about certain ideals preserved under dynamic equivalence, that neither of 
the first two systems listed above are dynamically equivalent to any other control 
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system with J ^ K = Q. The proof of the final theorem gives exphcit dynamic 
equivalences between the last three systems above. 

Theorem 15. The control system corresponding to = with two control vari- 
ables is not dynamically equivalent to any other control system with J = K = to 
which it is not static equivalent. 



Proof. Suppose fl = gSQu}, where Q is given by (11.41, S is given by (9.2 1, G is 
a generic nonautonomous static equivalence, and iv is the following coframing for 

±3 = 0. 

UJ~^ ^ dt 
oji — dxi — ui dt 

^2 = dx2 — U2 dt 
W3 = dx3 

The coframing ft = SQuj would then look as follows. 
Ci-^ = dt 

^1 = {{GDI dxi + {GDI dui + du2) - {{GDI + (G°)? + U2) dt 

= {dx2 + {GDI dxi) - {{GDI ui + U2) dt 
nl = dx3 

Now notice that the algebraic ideal fj" is preserved by G. But all of our equiva- 
lences also preserve t, and hence dt. Therefore, if Aqo has the coframing 

n-^ = dt 

^1 = dyi -~ gi{y,v) dt 
^2 = dy2~ g2{y,v) dt 
^3 = dys- g3{y,v) dt 



we get that {fli, f)", ^^3} = {dyi,dy2,dy3} mod dt. Since this is an integrable ideal 
that contains fig = dxs, we can arrange through the appropriate choice of G that 
dy3 = dx3. Note that this automatically satisfies dCt^^ = mod ft 1,(12,^13 since 
dfl^ is identically zero. 

Therefore y3 — X3 = 0. What we have done is taken any strict dynamic equiv- 
alence to the system ±3 = with J — K — and altered it via static equivalence 
to a strict dynamic equivalence to itself. So any control system that is dynamically 
equivalent to is = with J = if = is in fact a dynamic equivalence to a system 
that is static equivalent to ±3 = 0. □ 

Theorem 16. The control system corresponding to X3 = 1 with two control vari- 
ables is not dynamically equivalent to any other control system with J = K = to 
which it is not static equivalent. 

Proof. The proof is nearly identical to that of the previous theorem. Replace X3 = 
with X3 — 1, and proceed in the same fashion. □ 
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Note that the method used in the previous two theorems could also be applied 
to the case of ^3 = a;2. A difference occurs, however, when reaching the step 
ys =3^3 = X2- Since X2 is not necessarily equal to 2/2, we see that the resulting 
system may or may not necessarily be static equivalent to the original system 
xz = X2- It in fact turns out. as stated in the next theorem, that this new system 
need not be static equivalent to the original system. 

Theorem 17. The control systems is = X2, X2U1, l+a;2Ui are strictly dynamically 
equivalent to each other. 



Proof. The following sets of maps between infinite jet bundles give explicit dy- 
namic equivalences for the three systems. We will demonstrate that the maps take 
solutions of one control system to solutions of the other. The fact that the maps 
composed with their respective inverses are in fact the identity on solutions is simple 
enough and is left to the reader. 



Xi = Ui 


yi 


= Vi 


Zl = 




X2 = U2 


yi 


= V2 


Z2 = 


W2 


±3 = X2U1 


ya 


= 2/2 


23 = H 


h -22^1 



Equivalence maps: (x, u) -f^ (y, v) 

ip{x,u,u,...) 
^p~'^{y,v,v,...) 



= ( a;ia;2 - X3, U2, X2, X1U2, U2, ■■■) 

= {vi/y2, ys, y3Vi/y2-yi, 5 ' 2/2, ••• ) 

2/2 



Verifying solutions: 



yi = HX1X2-X3) ^1 ~ dt\y2 

= X1U2 _ y2Vi - V1V2 

= vi ~ y2^ 

= Ui 

X2 = m 



2/2 

V2 



2/2 = U2 
= V2 



h = X2 dt V 2/2 

= U2 , y2Vl - V1V2 

= Vl+2/3 o Vi 

= 2/2 2/2^ 



Equivalence map: (z, w) -H- (y, v) 



4>{Z,'W,U,...) = { Z3 ~ Z1Z2, W2, Z2, 1 - Z1W2, W2, ■ ■ ■ ) 
I -I, . s , 1-Vl Ul 1^2 - f2 - 2/21^1 

{y,v,w,...) = ( ,2/3,2/1+2/3 , 2 '2/2, 

2/2 2/2 2/2 
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Verifying solutions: 
d 



yi 



y2 



2/3 



dt - ''''^ 

1 + Z2'Wi — W\Z2 — ZiW2 



= W2 

= V2 

= Z2 

= W2 

= y2 



Zi = 



Z2 = 



Equivalence map: (x, u) -n- (z, w) 

9{X,U,U, ...) = ( Xi, X2, — 

U2 U2 

0~'^{Z,W,W,...) = { ^ - Zi, Z2, 



W2 



W2 



Note that 6* = V'"^ o ^. 
Verifying solutions: 



xi = 



X2 



1 

dt \'W2 
— W2 

Ul 
Z2 
W2 
U2 

d^ I' Z2_ 
dt \W2 



- Zi 



Wi 



Z3 



zz. 



Z2 y-wi 

X2U1 



W2 



d 1 — vi 

dt 1/2 

V1V2 -V2- y2Vl 



2/2^ 



y2 

W2 

d 
dt 



^3 = 37 h/i + 2/3 



2/2 



1 + 2/3 



{Vi - 1)V2 - 2/2*1 



= 1 + 2:2^1 



2/2^ 



- 0:3, -Ul - 



U2_ 
2 ' 



W2 



U2, ... ) 



- Z3, -Wi 2' ^2, ... ) 



Wo 



Z\ = 



Z2 



1 

dt \U2 
-U2 

^2^ 
Wi 

X2 
U2 

W2 

d /X2 
dt \U2 

1+X2 



1 + Z2WI 



- Xi 



X3 



-Ul 



U2_ 
W2^ 



□ 



13. Conclusions 

Wilkens showed that there are five equivalence classes of affine linear control 

systems with three state variables and two control variables under static equiva- 
lence. Below is a listing of how these classes combine using dynamic equivalence 
through one prolongation, using the static class representations presented in Elkin. 
Each equivalence class under dynamic equivalence is numbered. These nontrivial 
equivalences (or non-equivalences) are the work of this paper. 
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1 


Xl — Ui 
±2 = U2 

X3 = X2 


Xl = Ml 
X2 = U2 
±3 = X2U1 


Xl = Ml 
X2 = M2 
±3 = 1+ X2U1 


2 


Xl = ui 

±2 = U2 
±3 = 






3 


Xl = ui 

±2 = U2 

X3 = 1 







Future avenues of research into the classification of affine linear control systems 
under dynamic equivalence include looking at higher order equivalences (J and/or 
K > 0) as well as increasing the number of state and control variables. One 
obstacle to overcome with higher order equivalences and larger numbers of variables 
is that, unlike the case presented here where a unique <S exists, the problem quickly 
splits into many cases with different S. In addition, this method relies on the 
fact that affine linear systems in this dimension have already been classified under 
static equivalence, and the static equivalence problem for afiine control systems 
has only been completed in a few low-dimensional cases. Nevertheless, the further 
exploration of this decomposition may still yield new insights into the phenomenon 
of dynamic equivalence in general. 
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